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Abstract 

We reconsider the role that the possible detection of lepton flavour violating (LFV) decays of 
supersymmetric particles at the Large Hadron Collider (LHC) can play in helping reconstruct the 
underlying neutrino mass generation mechanism within the simplest high-scale minimal supergrav- 
ity (mSUGRA) seesaw schemes. We study in detail the LFV scalar tau decays at the LHC, assuming 
that the observed neutrino masses arise either through the pure type-I or the simpler type-H seesaw 
mechanism. We perform a scan over the mSUGRA parameter space in order to identify regions 
where lepton flavour violating decays of X2 can be maximized, while respecting current low-energy 
constraints, such as those coming from the bounds on Br(/i 67). We estimate the cross section 
for ~^ Xi+T + IJ'- Though insufficient for a full reconstruction of the seesaw, the search for LFV 
decays of supersymmetric states at the LHC brings complementary information to that coming 
from low energy neutrino oscillation experiments and LFV searches. 
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I. INTRODUCTION 



Neutrino oscillation experiments [l| have provided the first signal of physics beyond the 
standard model (SM). These measurements show that (a) neutrinos have a non-zero mass 
and (b) lepton flavour is violated. So far there is no experimental data that indicates that 
lepton number is also broken. However, one expects neutrinos to be Majorana particles, 
their mass at low energy being described by a unique (AL = 2) dimension-5 operator ^ 



Cm,. = ^{HL){HL). 



(1) 



where A is some high energy scale, / a dimensionless coupling constant, and H (L) the 
Higgs boson (lepton) doublets, respectively. Many model realizations of eq. JT]) exist, the 
most famous being the seesaw mechanism. The latter can be iinplemented via the exchange 
of a heavy singlet fermion, usually called type-I seesaw The exchange of a scalar 



triplet I lUj is now known as type-II seesaw ^ The exchange of a fermionic 

triplet [l2| is called type-Ill seesaw in [l3|. A list of generic 1-loop realizations of eq. ^ can 



also be found in 
are discussed in 



13| . Further seesaw realizations, such as the inverse and the linear seesaw, 
3. 

At "low" energies one cannot decide whether tree-level or loop physics generates eq. ([T]), 
nor can any measurements of neutrino angles, phases or masses distinguish between the dif- 
ferent tree-level seesaw realizations. Under the assumption of a pure type-I or pure type-II 
minimal supergravity seesaw mechanisms, we reconsider here the prospects for reconstruct- 
ing the underlying high energy parameters from a combination of different measurements. 
Clearly, observables outside the neutrino sector are needed in order to ultimately learn about 
the high energy parameters characterizing the seesaw. The classical tree-level realizations 
of the simplest type-I seesaw mechanism, unfortunately, can not be put to the test in a 
direct way. This can be straightforwardly understood by inverting eq. ([I]), which results in 

If the CERN LHC, due to take first data, finds signs of electroweak scale supersymmetry, 
indirect insight into the high-energy world might become possible through the search for 



fiavour violation effects [l5|,ll6|. Starting from fiavour diagonal soft supersymmetry breaking 
terms at some high energy "unification" scale, fiavour violation appears at lower energies due 
to the renormalization group evolution of the soft breaking parameters If the seesaw 

mechanism is responsible for the observed neutrino masses, the neutrino Yukawa couplings 
leave their imprint in the slepton mass matrices as first shown in jl^- Potentially large 
LFV is then induced by the fiavour off-diagonal structure in the Yukawa couplings required 
by the large mixing angles observed in oscillation experiments [20|. Expectations for LFV 



decavs such as L_— + 7 and k —>■ 31 j in the supersymmetric seesaw have been studied in 
261 . For the related process of /i — e conversion in nuclei see, for example 



21 
27 




The potential of LHC experiments in probing the allowed seesaw parameters 



^ This is the opposite convention to that used in 
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through measurements of masses and branching ratios of supersymmetric particles has been 
also discussed in Refs. 



29, 



31 



31 



32, 



33 



In two previous studies [3J, we have pointed out that ratios of branching ratios are 



especially useful for learning about the unknown seesaw parameters. In 
type-I seesaw was discussed, whereas 



35 



34 the case of 



addresses the case of seesaw type-II. For the type- 
I seesaw, there are in general too many unknown parameters that preclude making any 
definite predictions for LFV decays. In contrast, in the simplest type-II seesaw model (with 
only one triplet coupling to standard model leptons) neutrino mixing angles can be related 
to ratios such as Br(f2 — > eXi)/Br(f2 /^Xi)- 

It has been shown that, to a good approximation, such ratios do not depend on the 
mSUGRA parameter values. However, from an experimental point of view, calculations of 
absolute event rates are needed, before ratios of different final state channels can be studied. 
In |34l . |35| we took as reference just a few benchmark mSUGRA points, for which we have 



made detailed studies. In the present paper we calculate branching ratios and event rates 
over a large region of mSUGRA parameter space, in order to identify the maximal number 
of events one can expect in experiments at the LHC, while still respecting all low-energy 
constraints. 

The rest of the paper is organized as follows. In the next section we give a short summary 
of the theoretical setup. In section III we describe our numerical procedure and present our 
results. Finally we close in section IV with a discussion and a short summary. 



II. THEORY SETUP 

In order to fix the notation, we will briefiy discuss the main features of the seesaw mech- 
anism and mSUGRA. The type-I supersymmetric seesaw consists in extending the particle 
content of the MSSM by three gauge singlet "right-handed" neutrino superfields. The lep- 
tonic part of the superpotential is then 

W = Yj%HdE] + Yl%H^N'r + M.iV^iV;, z, j = 1, . . . , 3, (2) 

where and Y^, denote the charged lepton and neutrino Yukawa couplings, while Nf are the 
"right-handed" neutrino superfields with Mj Majorana mass terms. One can always choose 
a basis in which the Majorana mass matrix of the "right-handed" neutrinos is brought to 
diagonal form Mr = diag(Mi, M2, M3). Without loss of generality we will also assume that 
eq. ([2]) is written in the basis where the charged lepton Yukawa matrix is already diagonal. 
In this simple setup, the type-I seesaw model, as defined by eq. ([2]), is characterized by a 
total of 21 parameters, from which only 12 are measurable in the low-energy theory, as we 
discuss below. 

The effective mass matrix of the "left-handed" neutrinos at low energies is then given as 

= -^Yj . ■ r„ (3) 
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we 



so that, for each "right-handed" neutrino, there is one non-zero eigenvalue in m^. In eq. 
use the notation (Hu^d) = for the vacuum expectation values of the neutral components 
of the Higgs boson doublets. 

The parameters of eq.([2|) are defined at the Grand Unified Theory (GUT) scale, whereas 
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the entries of eq. (|3]) are measured at low energies. In order to connect t 
numerically solve the full set of renormalization group equations (RGE) 

Being complex symmetric, the light Majorana neutrino mass matrix in eq. 
nalized by a unitary 3x3 matrix U 



lese two scales we 



m,. 



is diago- 



(4) 



Inverting the seesaw equation, eq. ([3]), allows to express Y^, as [i 



n = V2- 



Vu 



Mn-R- 



(5) 



where is the diagonal matrix with rrii eigenvalues and in general i? is a complex orthogonal 
matrix. Note that, in the special case R = 1, Y^, contains only "diagonal" products y/Miirii. 
In this simplified case the 18 parameters in Yi, are reduced to 12. Note that in general type-I 
seesaw schemes, the unitary matrix diagonalizing the effective neutrino mass matrix differs 
from the lepton mixing matrix by terms of order D/Mji, where the D = Y^Vu- For the 
high-scale schemes considered here one can safely neglect these deviations ^. In this case 
we can set the diagonalization matrix as the lepton mixing matrix (partially) determined in 
neutrino oscillation measurements. 

Implementing the type-II seesaw mechanism within supersymmetry requires at least two 
SU{2) triplet states T12. A scalar triplet with mass below the GUT scale changes the 
running of gi and g2 in an unwanted way and gauge coupling unification is lost. If instead 
one adds only complete SU{b) multiplets (or GUT multiplets which can be decomposed 
into complete SU{b) multiplets) to the standard model particle content, the scale where 
couplings unify remains the same (at one loop level), only the value of the GUT coupling 
itself changes 
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Our numerical calculation uses an SU{h) inspired model [39|, which adds a pair of 15 
and 15 to the Minimal Supersymmetric Standard Model (MSSM) particle spectrum. This 
variant of the type-II seesaw mechanism, as discussed above, allows us to maintain gauge 
coupling unification even for -C M^, Mq being the unification scale. Under SU{?>) x 
SUl{2) X Uy{1) the 15 decomposes as 



15 

S 



S + T + Z 
(6,1,-^), 



(6) 



r~ (1,3,1), 



z 



(3,2.i: 



^ However for other type-I schemes, like the inverse seesaw [26|, |28| this approximation fails and leads to 
large LFV from right-handed neutrino exchange, even in the absence of supersymmetric contributions. 
For a systematic perturbative seesaw diagonalization method that covers all cases see the second paper 
in Ref. 01. 



4 



T has the correct quantum numbers to generate the dimension-5 operator of eq. ([T]). The 
SU (5) invariant superpotential reads 

W = -^Y^/ 5, ■ 15 ■ 5j + ^Ai5h ■ 15 ■ 5h + ^^25// ■ T5 ■ 5^ + Yf^ lOi ■ 5j ■ 5h 
+ Y^o lOi ■ lOj ■ 5h + A'/i5l5 ■ 15 + M^Eh ■ 5^ . (7) 

Here, 5 = {d^,L), 10 = (u^e^Q), 5^ = (t, ifa) and 5j^ = {i,Hi). Below the GUT scale, in 
the SU (5)-broken phase, the superpotential contains the terms 

i=(l^ L,T,L, + If drSd^) + If d^ZL, + d1Q,H, + Y^^ u^Q.H^ + 1^^ elL.H, 

+^{\^HdT^Ha + \2HuT2Hu) + MtTiT2 + M^ZiZa + MsS^S^ + . (8) 

As long as Mz ~ M5 ~ M^, gauge coupling unification will be preserved. Note that exact 
equality is not required for a successful unification. In our numerical studies we have taken 
into account the different running of these mass parameters. 

Integrating out the heavy triplets at their mass scale, the dimension-5 operator of eq. ^ 
is generated and after electroweak symmetry breaking the resulting neutrino mass matrix 
can be written as 

As in the case of the type-I seesaw, eq. ^ depends on the energy scale. In order to compute 
the neutrino mass rriy measured at low energies, one needs to know A2, Yt and Mt as input 
parameters at the high energy scale. As will be discussed in section IIIIl one can use an 
iterative procedure in order to find the high scale parameters from the low energy measured 
quantities. 

Note that, without loss of generality, we have the freedom to write eqs. (Ej) and ([7]) in the 
basis where the charged lepton mass matrix is diagonal, fitting the corresponding Yukawa 
couplings so as to reproduce the three measured charged lepton masses. However note that 
there are important differences between the type-I and type-II seesaw schemes. For example, 
in contrast to type-I, in a pure type-II seesaw scheme the unitary matrix U that diagonalizes 
eq. (jni) coincides with the lepton mixing matrix studied in neutrino oscillations. Moreover, in 
sequential type-I seesaw for each "right-handed" neutrino added there is one non-zero light 
neutrino mass eigenstate ^ . In contrast, in type-II seesaw one can produce three neutrino 
masses with just one pair of triplet superfields, with only one triplet directly coupling to 
leptons. This implies that in the minimal type-II seesaw one has less parameters than in 
the sequential type-I seesaw. Indeed from the 12 parameters in the complex symmetric Yt 



^ We do not consider here the possibility of having just two right-handed neutrino states in the type-I seesaw, 
called (3,2) in Ref. 0|. This could well account for the current neutrino data with just 12 parameters, 
instead of the 18 characterizing the sequential (3,3) seesaw considered here. 
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matrix, one can remove 3 phases by redefining the charged leptons [6|. In addition, from the 
3 complex parameters Ai,2 and M^, one does not enter, as only one of the triplets couples 
to leptons, and finally, two of the three phases can also be removed by field redefinitions. 
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The net result is that there are only 11 physical parameters governing neutrino physics 
This number is substantially smaller than the 18 free parameters describing the simplest 
type-I seesaw scheme containing three "right-handed" neutrinos ji^ 

At low energies a maximum of 9 neutrino parameters can be fixed by measuring lepton 
properties: 3 neutrino masses, 3 mixing angles and 3 CP phases. Thus from neutrino data 
only, neither type-I nor type-II seesaw schemes can be completely reconstructed, even in 
their simplest realizations. However, especially important in the following is the fact that 
low-energy neutrino angles are directly related to the high-energy Yukawa matrix in the 
type-II seesaw, whereas no such simple connection exists in the seesaw type-I (see also the 



discussion in 



|4H). 

As already commented, to a good approximation the lepton mixing matrix may be taken 
in unitary form, with three mixing angles and three physical CP phases (pij j^. Of 
these only the leptonic analogue of the Kobayashi-Maskawa phase 5, taken as the invariant 
combination 5 = 0i2 — (pn + 023 would enter the class of LFV processes discussed in this 
paper, so that we get the standard form. 



^ C12C13 S12C13 Si3e \ 

U = — S12C23 — Ci2S23Si3e* C12C23 — Si2S23Si3e* S23C13 

\ S12S23 - Ci2C23Si3e*'^ -C12S23 - Si2C23Si3e*'' C23C13 J 



(10) 



where Sij = sin 6ij, Cy = cos 6ij. Since no current experiment is sensitive enough to probe 
leptonic CP violation we take, for simplicity, 6 = 0. Neutrino oscillation experiments can 
be fitted with either a normal hierarchical spectrum (NH), or with an inverted hierarchy 
(IH) one. If one does not insist in ordering the neutrino mass eigenstates m,^., i = 1,2,3 
with respect to increasing mass, the matrix U can describe both possibilities without re- 
ordering of angles. In this convention, which we will use in the following, m^^^ ^ (m^g ~ 0) 
corresponds to normal (inverse) hierarchy and Su, ■S13 and S23 are the angles in both types 
of spectra. Basically su is measured in solar + reactor experiments, S23 in atmospheric + 
accelerator experiments and S13 is constrained by reactor neutrino oscillation data. 

In the general MSSM, LFV off-diagonal entries in the slepton mass matrices involve 
additional free parameters which arise from the mechanism of supersymmetry breaking. In 
order to relate LFV in the slepton sector with the LFV encoded in Y^, or Yr one must 
assume some particular scheme for supersymmetry breaking. For simplicity and definiteness 
we will adopt mSUGRA boundary conditions, characterized by four continuous real and one 
discrete free parameter, usually denoted as 

mo, Ml/2, Ao, tan/3, Sgn(/i) . (11) 



^ We are treating the three charged lepton masses as experimentally determined parameters. 
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Here, mo is the common scalar mass, M1/2 the gaugino mass and Aq the common trihnear 
parameter, all defined at the grand unification scale, Mq — 2 ■ 10^^ GeV. The remaining two 
parameters are tan/? = v^/vd and the sign of the Higgs mixing parameter fi. 

In order to have a qualitative understanding of the magnitudes of the LFV rates we first 
present approximate leading-log analytical solutions for the renormalization group equa- 
tions ^. For the case of type-I seesaw, the LFV elements induced in the charged left-slepton 
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mass matrix by renormalization group evolution can be approximated as 

(AM?),, = -^(3m2 + , (12) 

where Yiy is given in terms of the neutrino parameters by eq. ([5]) and the factor L is defined 



as 

^fe; = log(^)4z. (13) 

Similarly, one can get an analogous approximate expression for the off-diagonal eleinents 
of the charged left-slepton mass matrix characterizing LFV in type-II seesaw schemes 



39 . 



III. NUMERICAL RESULTS 

Due to the non-trivial structure of the neutrino Yukawa matrix Y^, in eq. ([5]) and of 
Yt in eq. (I9|) for type-I and type-II seesaw, respectively, the slepton mass matrices contain 
calculable LFV entries In order to determine their magnitude we solve the complete 



set of renormalization group equations, given in [22|, l36|, l39|| . All results presented below have 
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been obtained with the lepton fiavour violating version of the program package SPheno 
where the RGBs for the MSSM part are implemented at the 2-loop level. For definiteness we 
set neutrino mass squared differences to their current best fit values [20| and fix the angles 
to the Tri-Bi-Maximal (TBM) values H. 

Fixing the values of other mSUGRA parameters, we used SPheno to perform a numerical 
scan over the mo-Mi/2 plane. For each point in this plane, we adjust the value of Mr (Mr) 
in order to keep the low energy LFV observable BR(/x — > 67) within its present exp erimental 
upper bound or within the expected sensitivity of the upcoming experiments 



For type-I seesaw our numerical procedure to fit these masses is as follows. As we have 
already commented, the large number of free parameters characterizing even the simplest 
type-I seesaw schemes forces us to make simplifying assumptions in inverting the seesaw 
equation, eq. As a first step we assume degenerate "right-handed" neutrinos and the 
simplest possible, flavour-blind, structure for the matrix R, i.e. 

i? = 1, Mmj = Mr 6,, . (14) 



^ Note that in the numerical code that leads to the results presented in our plots we have numerically solved 
the full set of RGEs. 
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Moreover, we fix the values of the light neutrino masses and Yukawa couplings to reproduce 
the TBM angle values. In order to determine the resulting LFV observables we numerically 
integrate the renormalization group equations taking into account the flavour structure of 
the Yi, matrix. We integrate out every "right-handed" neutrino and its superpartner at the 
scale associated to its mass, and calculate the corresponding contribution to the dimension- 
five operator which is evolved to the electroweak scale. This way we obtain the exact 
neutrino masses and mixing angles for this first guess. The difference between the results 
numerically obtained and the input numbers is then minimized in an iterative procedure 
until convergence is achieved. 

For the type-II seesaw the calculations are performed for the 15-plet case, under the 
assumption Yz = Yt = Ys at Mq, as discussed above, and including the one-loop RGEs 
for the new parameters in SPheno. For consistency, we have also included 1-loop threshold 
corrections for gauge couplings and gaugino mass parameters at the scale corresponding to 
the mass of the triplet, M^. The MSSM part is implemented at the 2-loop level and, thus, 
in principle one should also consistently include the effect of the 15-plets for all parameters 
at this level. However, as discussed in the correct fit of the neutrino data requires 

that either the triplet (15-plet) Yukawa couplings are small and/or that Mt is close to Mq, 
implying that the ratio Mt/Mq is significantly smaller than Mc/rriz and thus one expects 
only small effects. Inverting the seesaw equation for any fixed value of A2 in eq. ([9]), one can 
get a first guess of the Yukawa couplings for any fixed values of the light neutrino masses 
as a function of the corresponding triplet mass. This first guess will not give the correct 
Yukawa couplings, since the neutrino masses and mixing angles are measured at low energy, 
whereas for the calculation of we need to insert the parameters at the high-energy scale. 
However, we can use this first guess to numerically run the RGEs to obtain the exact neutrino 
parameters (at low energies) for these input values. The difference between the results 
obtained numerically and the input can then be minimized in a simple iterative procedure 
until convergence is achieved. As long as neutrino Yukawas are not large, convergence is 
reached in a few steps. However, in type-II seesaw schemes, the Yukawa couplings run 
stronger than in the type-I seesaw, thus our initial guess can sizeable deviate from the exact 
Yukawa coupling values. Since neutrino oscillation data requires at least one neutrino mass 
to be larger than about 0.05 eV, we do not find any solutions for Mt > 10^^ GeV. 

Finally, the calculation of cross sections for the production of supersymmetric particles 
was done using Prospino j^l- The input data was taken from SPheno using the SUSY Les 
Houches Accord standard format (4^ . 

A. Br(/i ^ e + 7) for type-I and type-II 

In Fig. [Hand Fig. [2] we show the contours for Br(;U — > e + 7) in the ttiq, M1/2 plane for pure 
type-I and type-II seesaw schemes, respectively. On the left panel of Fig. [T] we chose a low 
value for Mr = 10^^, while on the right panel a value of 10^^ GeV was chosen. In Fig. [2]the 
same dependence is shown for the type-II seesaw mass scale Mt- This shows the dependence 
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FIG. 1: Contours of Br(/x ^ e + 7) in the mo,Mi/2 plane for our standard choice of parameters: 
/i > 0, tan/3 = 10 and Aq = GeV, for type-I with degenerate RH neutrinos. On the left panel 
Mr = 10^3 QgY^ on the right panel Mr = 10^'' GeV. 




FIG. 2: Contours of Br(// ^ e + 7) in the mQ,Mi/2 plane for Ai = 0.02 and A2 = 0.5 and for our 
standard choice of parameters: > 0, tan/3 = 10 and Aq = GeV, for type-II seesaw. On the left 
panel Mr = 10^^ GeV, on the right panel Mr = 10^^ GeV. 



of LFV rates on the choice of scale Mr (Mr). The complicated features displayed on these 
plots are due to cancellations between the chargino and neutralino amplitudes contributing 



to — > e + 7, as is well known [23|, |3j|. First we note that large parts of parameter space 



fall within the regions of sensitivity of upcoming experiments like MEG 



44 . The contours 
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for BR(/i —^67) are deformed in type-II as compared with respect to those for type-I 
seesaw. The reason for this is that the addition of gauge non-singlet states in type-II seesaw 
increases the dependence on the renormalization scale of the neutrino Yukawa coupling and 
also affects the supersymmetric spectrum, which alters the region where BR(/i — > 67) is 
strongly suppressed. For our subsequent discussion, the most important point is that for 
each point in the mo,Mi/2 plane there will be a maximum value of Mji (Mt) that will 
give the maximum possible rates of LFV compatible with current experimental bounds, 
BR(/i ^ 67) < 1.2 X 10-" 
experiments like MEG, BR(/i 



42|, or with expected sensitivities to be reached at upcoming 



6 7) < 10-13 [44| 



B. LFV stau decays 



The eagerly awaited production of supersymmetric particles at the LHC would open new 
opportunities for the study of flavour violation in the supersymmetric sector Here we 
study how the LFV decays of staus may provide valuable cross-checks of neutrino properties 
determined at low energies as well as complementary information on the origin of neutrino 
mass. 




FIG. 3: Br(f2 fJ- + Xi) (left panel) and Br(f2 — > e + Xi) (right panel), in the mo,Mi/2 plane for 
our standard choice of parameters: > 0, tan/? = 10 and Aq = GeV, for type-I seesaw, imposing 
Br(// ^ e + 7) < 1.2 • 10-1^ 

The expected LFV branching ratios for T2 — > + Xi ^ind f2 e + Xi depend on the 
choice of the mSUGRA parameters. After a full scan over the mSUGRA parameter space 
we found that the dependence on Aq and on the sign of fi is weaker, but that the rates 
decreased with increasing values of tan p. Therefore, we chose our standard point with a 
relatively low value of tan/? = 10, and for deflniteness took > 0, and Aq = 0. In Fig. [3] 
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FIG. 4: Br(f2 ^ /i + X?) (left panel) and Br(f2 
our standard choice of parameters: > 0, tan (3 = 10 and Aq 
Br(// ^ e + 7) < IQ-^^^ 



(right panel), in the mQ,Mi/2 plane for 
GeV, for type-I seesaw, imposing 



we show the contour plots for the LFV decays T2 — > /i + Xi (left panel) and T2 — ^ e + x? 
(right panel) in the mo,Mi/2 plane for our standard choice of mSUGRA parameters for 
the simplest pure type-I seesaw scheme. One sees that there are regions in parameter 
space where the LFV decays of the f2 can be as large as of order 10~^. In these plots the 
values of Mr were chosen as to obtain the maximum LFV compatible with the present 



experimental limit of Br(/i — > e + 7) < 1.2 ■ 10~^^ [42|. Also shown in these plots are the 
exclusion regions coming from the LEP constraints on SUSY masses and also the exclusion 
obtained when the neutralino is not the LSP ^. In Fig. [4] we show the same contour plots 
for Br(/i — e + 7) < 10^^'^, which will be achievable in the coming experiments 4j|. Also 
in this case one observes in Fig. [4] that the LFV stau decay rates may exceed the 10% level. 
Notice also that the nontrivial features present in in Figs. [3] and Fig. [4] reflect the well-known 
cancellations between chargino and neutralino contributions to /i — > e + 7 already discussed 
above. 

In Fig. [5] and Fig. [6] the same type of plots are shown for type-II seesaw. A comparison 
of these figures shows that, qualitatively, the behavior is very similar for the two types of 
seesaw. In both cases, the larger rates for r2 ^ e + Xi are more constrained in parameter 
space than those for f2 + Xi- Notice however that there is an important difference 
between type-I and type-II seesaw, coming from the presence of the Higgs triplets that 
contribute sizeably to the running of the type-II beta functions. This gets reflected in the 
supersymmetric particle spectra and hence in the shapes of the red (shaded) regions in 



Note that we did not display the constraints coming from Dark Matter (DM) relic abundance. 
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Fig. [5] and Fig. [6l One can observe, indeed, that the regions where the stau is the lightest 
supersymmetric particle, as well as the regions already excluded by LEP2 are substantially 
different for type-II seesaw, as compared to the corresponding ones for type-I. This follows 
from the modification in the beta functions introduced by the addition of the Higgs triplets, 
making Mi and M2 smaller in type-II than in type-I seesaw for the same value of M1/2. The 
variation with the mSUGRA parameters is illustrated in Fig. [7] (type-I) and Fig. [8] (type-II) 
for the parameter Aq and in Fig. [9] (type-I) and Fig. [10] (type-II) for tan /5. We can see that 
there is not much variation with Aq, while the rates decrease rapidly with increasing values 
of tan/3. The reason for this is that BR(/i — > e + 7) increases as tan^/5, thus constraining 
more strongly the maximum attainable stau LFV rates. This effect is stronger for type-I 
as can be seen by noting the different values for the contour levels in Fig. [9] and Fig. [TOl 
The variation with the sign of /i is weak and we do not show it here. So, in summary, large 
LFV rates prefer moderate values of tan j3 and this explains a posteriori the choice of our 
standard parameters. 




FIG. 5: Br(f2 ^ /i + Xi) (left panel) and Br(f2 — > e + x?) (right panel), in the uiq, M1/2 plane for 
Ai = 0.02 and A2 = 0.5 and our standard choice of parameters: /i > 0, tan /3 = 10 and Aq = GeV, 
for type-II seesaw, imposing Br(/i — > e + 7) < 1.2 • 10^^^. 



C. Total production cross section of X2 

As important as having a large branching ratio into a LFV final state, is to be able 
to produce a large enough event sample. In order to estimate the number of LFV events 
expected at the LHC, one notes that, from Figs. [3] - [IDl in the regions where the LFV is 
sizeable, the direct production of staus at the LHC is negligible compared to that which arises 
from cascade decays of heavier neutralinos, mainly X2- We focus on the X25 because decays 
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FIG. 6: Br(f2 ^ fi + x?) (left panel) and Br(f2 ^ e + Xi) (right panel), in the tuq, M1/2 plane, for 
Ai = 0.02 and A2 = 0.5 and our standard choice of parameters: /i > 0, tan/3 = 10 and Aq = GeV, 
for type-II seesaw, imposing Br(/i — > e + 7) < 10"^'^. 




FIG. 7: Br(f2 — > /i + Xi) (left panel) and Br(f2 — > e + x?) (right panel), in the mo, M1/2 plane for 
standard choice of parameters: /U > 0, tan/3 = 10 but different Aq = —300 GeV, for type-I seesaw, 
imposing Br(// — > e + 7) < 1.2 • 10^^^. 



such as X2 ~^ f^TXi ctre sensitive to flavour violation, whereas in the corresponding chargino 
decays the flavour information is lost. Hence we flrst compute the total X2 production cross 
section. In the left panel of Fig. [H] we show the results for the cross section for X2 production 
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FIG. 8: Br(f2 + X?) (left panel) and Br(f2 e + x?) (right panel), in the mQ,Mi/2 plane 
for Ai = 0.02 and A2 = 0.5 and standard choice of parameters: /x > 0, tan/3 = 10 but different 
Aq = —300 GeV, for type-II seesaw, imposing Br(/i — > e + 7) < 1.2 • 10"^^. 



BR(t: ^ |i x°), tan|3=30, A(,=0 (GeV) 
BR(n^ e y)=1.2 X 10"''^ ; |i > 
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mo (GeV) 

FIG. 9: Br(f2 ^ /i + X?) (left panel) and Br(f2 
standard choice of parameters: fi > 0, Aq = b 
Br(// ^ e + 7) < 1.2 • 10-1^ 



BR(x ^ e x°), tan|3=30, Ao=0 (GeV) 
BR(n^ey)=1.2x 10"''^ m>0 




100 200 300 400 500 600 700 800 900 
mo (GeV) 

e + Xi) (right panel), in the mo, M1/2 plane for 
different tan/3 = 30, for type-I seesaw, imposing 



as a function of M1/2, for different choices of mo and for our standard choice of mSUGRA 
parameters: /i > 0, tan/3 = 10 and Aq = GeV, for the pure type-I mSUGRA seesaw 
scheme. This choice of mSUGRA parameters corresponds, as will be discussed below, to 
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BR(t ^ |i tanp=30, Ao=0 (GeV) BR(t ^ e x°), tanp=30, Ao=0 (GeV) 

BR(n^ e y)=1.2x 10"''^ |i> BR(n ^ e y)=1 .2 x 10"''^ ; ^> 
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FIG. 10: Br(f2 ^ M + X?) (left panel) and Br(f2 ^ e + x?) (right panel), for Ai = 0.02 and As = 0.5, 
in the mo,Mi/2 plane for standard choice of parameters: fi > 0, Aq = 0, but different tan/? = 30, 
for type-II seesaw, imposing Br(/i — > e + 7) < 1.2 • 10^^^. 




[GeV] Mj^2 [GeV] 

FIG. 11: Production cross section (at leading order) of X2 versus for varying mo, and for our 

standard choice of parameters: /U > 0, tan /? = 10 and Aq = GeV, in type-I seesaw (left panel) 
and type-II seesaw (right panel) for Ai = 0.02 and A2 = 0.5. 



the case where the branching ratios of the LFV stau decays are the largest. This result 
was obtained using the Prospino code |46|| at Leading Order (LO) approximation. We have 
checked that the Next to Leading Order (NLO) calculation only changes the results slightly, 



due to an appropriate choice of the renormalization scale [46|. So, in all cross sections 
presented here, we only used the LO approximation. The corresponding results for type-II 
seesaw are shown in the right panel of Fig. [HI for the same choice of mSUGRA parameters 
and for Ai = 0.02 and A2 = 0.5. 
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D. Total production of X2 times BR to /i-r lepton pair 

In order to get an estimate of the expected number of LFV events at the LHC we now 
use a combination of the Prospino and SPheno codes to evaluate the product of the X2 
production cross section times the branching ratios into LFV processes. Once we know the 
luminosity at LHC we can multiply it with the above product to get the number of events. 

In Fig. [121 we have plotted, for type-I seesaw (left panel) and type-II (right panel), the 
production cross section at leading order of the second lightest neutralino cr{x2) times the 
BR of X2 goiiig to the opposite-sign dilepton signal x? /^t as a function of M1/2, for different 
values of mg. We have fixed the rest of the mSUGRA parameters to our standard mSUGRA 
point and imposed an upper limit on Br(yU — > e + 7) < 1.2 ■ 10~^^. In type-I seesaw, the 
number of events of the opposite-sign dilepton signal X2 ^ Xif^T can be of the order of 10^ 
for mo ~ 100 GeV and M1/2 ~ [450, 600] GeV, assuming a luminosity C = 100 fb~^. In 
type-II seesaw, there can be a maximum number of events of the order of 10^ for mo ~ 100 
GeV and M1/2 ~ [600, 800] GeV. 
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FIG. 12: Production cross section (at leading order) of X2 times BR of X2 going to fi-r lepton pair 
versus M1/2 for mo = 100 GeV (red), 200 GeV (green), 300 GeV (blue) and 500 GeV (magenta), and 
for our standard choice of parameters: /U > 0, tan /? = 10 and Aq = GeV, for type-I (left panel) and 
for type-II seesaw (right panel) with Ai = 0.02 and A2 = 0.5, imposing Br(/i — > e-l-7) < 1.2- 10~^^. 



For type-II seesaw where we have less parameters, we can look at variations of the re- 
sult with the values of the triplet Higgs boson coupling A2, a parameter that can not be 
determined from neutrino data alone as it appears only in the ratio A2/MT, see eq. Q. In 
Fig. Us] we show the dependence of the product of cross section times LFV branching ratios 
as function of A2 for our standard point. We should mention that the other Higgs boson 
triplet coupling Ai, does not contribute to LFV decays, and hence is left undetermined by 
this analysis. 



As has been discussed in |48||, the dominant standard model backgrounds for the process 
considered are expected to be WW and tt production. The cuts necessary to reduce this 
background will depend on the details of the SUSY spectrum and a detailed investigation is 
beyond the scope of this paper. The results of |48| suggest that the signal should be visible 



16 



_ 10' 

e 

^ lo' 
pa 

^ 10-' 

O IN 

^ 10"' 



10" 







— ?i,=0.5 














- ?i-=0.9 












[ \ 







































400 



M 



1/2 



600 
[GeV] 



800 



FIG. 13: Production cross section (at leading order) of X2 times BR of X2 goi^ig to /x-r lepton pair 
versus M1/2, for our standard choice of parameters: /i > 0, tan/? = 10 and Aq = GeV, for type-II 
seesaw, imposing Br(/i — > e + 7) < 1.2 • 10"^^, for a fixed value of mo = 100 GeV and different 
values of A2 = 0.1 (green), 0.5 (red), 0.9 (blue). 



for ct(x^)xBR of order C(10) fb. 



IV. CONCLUSIONS AND OUTLOOK 

Low energy neutrino experiments, including oscillation studies and neutrinoless double- 
beta decay searches may, optimistically, determine at most 9 neutrino parameters: the 
3 neutrino masses, the 3 mixing angles and potentially the 3 CP violating phases. This 
is insufficient to fully reconstruct the underlying mechanism of neutrino mass generation. 
Under the assumption that neutrino masses arise a la seesaw, we have considered the simplest 
pure type-I or pure type-II seesaw schemes in mSUGRA. 

We have performed a full scan over the mSUGRA parameter space in order to identify 
regions where LFV decays of X2 can be maximal, while still respecting low-energy constraints 
that follow from the upper bounds on Br(yU — > 67). We have also estimated the expected 
number of events for X2 ~^ Xi + T + f^: for a sample luminosity of £ = 100 fb~^. The expected 
number of events for the other channel X2 ~^ + + c is always smaller, as can be seen from 
the LVF branching ratios presented in section IIII B[ We have found that the pure seesaw- 
II scheme is substantially simpler and comes closer to be fully reconstructable, provided 
additional LFV decays are detected and some supersymmetric particles are discovered at 
the Large Hadron Collider. 

Note that in what concerns the expected maximum number of events both type-I and 
type-II schemes give similar results. However, as we have seen, given their smaller num- 
ber of parameters, type-II seesaw schemes are more likely to be reconstructable through a 
combination of low energy neutrino measurements, with the possible detection of supersym- 
metric states and lepton flavour violation at the LHC. This should encourage one to perform 
full-fledged dedicated simulations, in order to ascertain their feasibility within realistic ex- 
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perimental conditions 
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Finally we note that we have not exploited the fact the LFV might induce new "edge 



49 . We have focused here on LHC, but mention 



variables", giving additional information 
that a future ILC would be much more suited for measuring LFV SUSY processes [29|, l50|, 
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